Thermal rectification in coupled rotor lattices with gradient mass 
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Abstract 

Heat conduction through one-dimensional (1D) coupled rotator lattices is investigated in the 
presence of mass gradient. It is found that thermal current in the direction of mass increasing is 
not asymmetric with heat flux though mass decreasing system, which is called thermal rectification. 
Moreover, we find that the larger is the mass gradient, the more evident is the thermal rectification 
effect, which is consistent with FPU lattice with mass gradient. Based on the influences of the 
thermal rectification on shape parameters K and A of nearest-neighbor interactions, the average 
temperature of atomic chain and the system size, a optimum thermal rectifier is designed, which 
thermal rectification efficiency has reached 70. These investigations would contribute to controlling 


and manipulating thermal current. 
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I. INTRODUCTION 


Heat conduction is an ancient physical question, which was described by Fourier Law as 
J = KAT where « is a constant, independent of system size. However, in recent years plenty 
of investigations have suggested that heat conduction in low dimensional systems diverges as 
the system size, for example FPU-a, FPU-{, diatomic Toda and so on{1]. These discoveries 
have attracted massive attention on heat conduction in low dimensional systems. Especial- 
ly, the proposal[2] of processing information by controlling and manipulating phonons(heat 
carrier) has motivated people to have designed various thermal devices, similar to electronic 
counterparts, thermal transistors|3], thermal logic gates|4], thermal memory[5] and thermal 
limiters|6]. To our delighted, the solid state thermal rectifier from carbon nanotubes(CNTS) 
with nonuniform axial mass distribution has been realized experimentally [7]. 

There exists various studies of designing various thermal rectifiers and negative differen- 
tial thermal resistance [8-25]. Li and his-coworker firstly presented the model|8] consisting 
of three segment FK chain of different interaction parameters to control the strength of the 
substrate potential. Later on, another models[9, 10] coupled two nonlinear one dimensional 
lattices increase the thermal rectification magnitude greatly. These investigations have been 
explained by mismatch/match of frequency spectrum at the two ends of the chain. Yang and 
his-coworker|19] theoretically showed thermal rectification in graded mass of FPU-8 lattices, 
and found that the larger is the mass gradient, the more obvious is the thermal rectification 
effect. Subsequently, Shah|20] studied the FPU-@ lattices with exponential mass graded 
lattices, and showed that the exponential mass graded material is a better and genuine 
choice for thermal rectification. Hu[21] explored the thermal rectifiacion in one dimension- 
al mass-graded harmonic lattice with harmonic and anharmonic on-site potential. It was 
found that thermal rectification occurs in the model with the anharmonic on-site potential 
and vanishes with the harmonic one, which suggests nonlinearity is necessary for thermal 
rectification. On the other hand, extensive investigations have done to determine the mech- 
anisms that lead to thermal rectification in graded mass systems. Nan and Wang found 
that thermal rectification results from the transmission of high-frequency phonons by ap- 
plying the nonequilibrium Green’s function method[24]. Emmanuel have showed the graded 
anharmonic systems as genuine thermal rectifier by applying self-consistent reservoirs|22]. 


Zhong and co-workers [25] found that graphene nanoribbons with thickness-asymmetry have 


a good thermal rectification and the thermal rectification factor depends on temperature as 
well as the thickness-ratio of the two-segment. 

In this paper, we investigate one-dimensional (1D) graded mass lattices with nonlinear 
nearest-neighbor interactions, and discover thermal rectification. In order to amplify the 
thermal rectification effect, the dependence of thermal rectification on the mass gradient, 
the shape parameter A and K, the average temperature To and the system size have been 


investigated. 


Il. MODEL AND METHODS 


The Hamilton of one-dimensional (1D) chain under consideration is of the form [1, 2] 


=F 


where x; is the instantaneous displacement of the ith particle from its equilibrium position, 


— a) + U(x), (1) 


P, denotes the ith particle’s instantaneous momentum, M; is the ith mass of the particle. 
V(axi41 — zi) is the nearest-neighbor interaction potential, U(x;) is the on-site potential. 
N is the total number of particles in the chain. In the coupled rotator model V takes an 


anharmonic form |1, 2] 
V = K(1 — cos A(z; — 2;_1)), U (z;) = 0. (2) 


For linear mass gradient, the mass of the ith particle is given by M; = Mmaz — (i — 1)(Mmar — 
Mmin)/(N — 1), where Maz is the mass of the particle at left end (first particle) and Min 
is the mass of the particle at right end (last particle) of the chain. In our simulation, 
Langevin heat baths[1] were imposed on the two ends and the fixed-boundary conditions 


Xo = £y+1 = 0 were adopted. Therefore the equations of motion are: 


A 
Miti = — — iti + Mi, 3 
y= Fo yit (3) 
where yi = (0:1 + ĉin) and Ni = Hi, + Nô: n. The equations of motion for the central 
particles (i = 2, 3..., N — 1) are, 


Miti = KA{sin A(ti41 = Ti) — sin A(x; a Li-1)]; (4) 


the equations of motion for i = 1 and 7 = N particles are given by, 


Mı =K Alsin A(z2 — £1) — sin A(zı — zo)| — 41 
+ nu(t), 
Mnin =K Alsin A(ry41 om zy) — sin A(xry = XLN_-1)| ) 
6 
where y is the friction coefficient. The dot stands for the derivation with respect to time t. 


The noise terms n+ (t) denote a Gaussian white noise that satisfies the fluctuation dissipation 


relations 
(nu(t)n(t)) = 2ykeTud(t — t), (7) 
(ni(t)ni(t )) = 2ykeT146(t — t'), (8) 


kpg denotes Boltzmann’s constant and we set kg=1 and y = 1 in our simulations. 

In this paper, we set Ty = To(1 + AT) and T; = To(1 — AT), where To = (Ty + Tr)/2 
is the average temperature and AT = Ty — Ty, is the temperature bias of the two ends of 
particle chain. We set Jz as the thermal current in the direction of mass decreasing when 
the high temperature bath (on the left) was imposed on Mmar (the maximum particle mass), 
and Jp as the heat flux in the direction of mass increasing when the high temperature bath 
(on the right) was coupled with Mmin (the minimum particle mass). 

For simplicity, the local heat flux is defined as j; = (tF (x; — xj_-1)), where F=-OV /Ox 
and the notation (...) denotes a steady-state average. After the system reaches a stationary 
state, J; is independent of site position i, so that the local flux can be denoted as 7 and the 
total heat flux J = Nj. In our simulations, the equations of motion(4-6) are integrated by 
using a second-order stochastic Runge-Kutta algorithm|26] with a small time step (h=0.002). 
The simulations are performed long enough to allow the system to reach a nonequilibrium 


steady state in which the local heat flux is a constant along the chain. 


Ill. NUMERICAL RESULTS AND DISCUSSION 


Figure 1 presents the dependence of the heat fluxes (Jr, Jz) (Jz denotes the thermal 
current from left to right and Jp is the heat flux from right to left) on the temperature 


gradient AT. It is observed that the heat fluxes increase as the temperature difference AT’ 
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FIG. 1: The heat flux J as a function of the temperature gradient AT. 
To = 0.5, A= 1.0, K = 1.0, Mmaz = 10.0, Mmin = 1.0, N = 32. 
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FIG. 2: (a)The heat fluxes (Jr, Jz) as a function of the mass gradient AM 


of the heat rectification efficiency R =| Jr/Jz | on the mass gradient AM. 


T = 0.5, A = 1.0, K = 1.0, Mmin = 1.0, N = 32. 


. (b) The dependence 


The parameters are 


increases. However, compared on Jz, the amplitude of Jr is much larger, which is called 


thermal rectification depicting thermal currents are not asymmetry when the heat baths are 


reversed. 


Figure 2(a) shows the dependence of the heat fluxes J on the mass gradient AM = 
Mimar/Mmin at Mmin = 1.0 for Ty = 0.5 and Ty; = 0.1. If the mass gradient AM is 


increased, heat currents (Jr, JL) through systems decrease monotonically. Obviously, the 


heat flux flows more easily in the direction of mass increasing for any values of mass gradient. 


It might be assumed: the uniform distribution of atomic masses and the nonlinear nearest- 
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FIG. 3: (a) The heat fluxes (Jg, Jz) as a function of the shape parameter A of the nearest- 
neighbor interaction. (b) the dependence of the heat rectification efficiency R =| JR/Jz | on the 
shape parameter A. The parameters are Tr = 0.1, Ty = 0.5, K = 1.0, Mmax = 10.0, Mmin = 1.0, 
N =32. 


neighbor interaction lead to the mismatch of power spectrum at the two ends of atomic chain. 
Therefore there exists asymmetric thermal currents. In order to describe quantitatively the 
thermal rectification efficiency, we introduce the ratio R =| JR/Jz |, when R = 1 shows no 
thermal rectification effect, and the larger is the value of R =| Jr/Jz |, the better is thermal 
rectification effect. 

Figure 2(b) shows the dependence of the rectification efficiency R on the mass gradient 
AM. As the mass gradient AM becomes larger, the thermal rectification efficiency R 
increases. It is found that the mass gradient AM contributes to thermal rectification, which 


is consistent with previous investigations[19—21]. 


Figure 3(a) depicts the dependence of the heat currents (Jr, Jz) on the shape parameter 
A. The heat fluxes (Jp, Jz)all increases with increment of the shape parameter A. It 
is observed easily that J; and Jg do not increase at the same pace, which suggests the 
shape parameter A could have effects on thermal rectification. Because the augment of the 
shape parameter A gives rise to a larger nonlinearity of the nearest-neighbor interaction. 
It improves the probability of the umklapp processes, which increases the degree of the 
mismatch of power spectrum of the two ends of atomic chain. Based on Figure 3(b), the 


more evident is the thermal rectification effect, the larger is the shape parameter A of 
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FIG. 4: (a)The heat fluxes (Jp, Jz) as a function of the strength K of the nearest-neighbor interac- 
tion. (b) the dependence of the thermal rectification efficiency R =| Jg/Jz | on the strength K of 
the nearest-neighbor interaction. The parameters are Tr = 0.1, Ty = 0.5, A = 1.0, Mmaz = 10.0, 
Mmin = 1.0, N = 32. 


nearest-neighbor interactions. 


Figure 4(a) shows the dependence of the heat fluxes J; and Jg on the strength K of 
the nearest-neighbor interaction at Ty = 0.5 and Tz = 0.1. Jz that is the heat flux in the 
direction of decreasing mass increases for K = [0.1, 1.0], and barely changes with increasing 
the strength K of the nearest-neighbor interaction. More interestingly, Jr that is the 
average thermal current in the direction of increasing mass decreases with the increment 
of the strength K of the nearest-neighbor interaction (further investigations see Fig. 5). 
Figure 4(b) describes that the thermal rectification ratio R as a function of the strength 
K of the nearest-neighbor interaction. It is found that the thermal rectification effect 
performs better at the minimum and maximum limit of the strength K of nearest-neighbor 


interactions. 


Figure 5(a) depicts the dependence of the thermal currents (Jp, J) on the average tem- 
perature To of chain. For a larger Tọ, both Jz and Jr decrease with increasing the average 
temperature To, which is related to nonlinear localized rotational modes of chain[27, 28]. 
Figure 5(b) presents thermal rectification efficiency R as a function of the average temper- 


ature Tọ. The rectification efficiency R monotonically increases as the average temperature 
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FIG. 5: (a)The heat fluxes (Jp, Jz) as a function of the average temperature To = (Ty + Tg)/2. 
(b) the dependence of the thermal rectification R =| JR/Jz | on the average temperature Ty. The 
parameters are AT=0.4, Ty = To(1 + AT), Tr = To(1 — AT), A = 1.0, K = 1.0, Mmax = 10.0, 
Minin = 1.0, N = 32. 
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FIG. 6: The heat fluxes J as a function of the strength K of the nearest-neighbor interaction for 
different values of the heaviest mass. (a) Mmaz=1.0, 3.0, 5.0 and 8.0. (b) Mmaz=10.0, 15.0, 20.0 
and 40.0. The parameters are A = 1.0, Ty = 0.1, Ty = 0.5, Mmin = 1.0, N = 32. 


To increases. It might be inferred: more phonons would be generated with rising the average 
temperature To. Large mounts of phonons of low-frequency collides and would be converted 
to phonons of high-frequency, which increases the weight of phonons of high-frequency and 


reduces the weight of phonons of low-frequency in entire phonons spectrum. 


In order to investigate the effects of the mass gradient on the phenomenon that Jp 
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FIG. 7: (a)The heat fluxes (Jr, JL) as a function of the total particles number N. (b) the depen- 
dence of the thermal rectification efficiency R on the system size N. The parameters are A = 10.0, 


K = 1.0, Mmaz = 10.0, Mmin = 1.0, Tr = 0.1, Ty = 0.5. 


decreases with increasing the strength K of nearest-neighbor interactions, we obtained Figure 
6(a) and Figure 6(b) by nonequilibrium molecular dynamics simulations(NEMD). Figure 6 
shows the dependence of Jr on the strength K of nearest-neighbor interactions for various 
values of mass gradient. When the mass gradient is small (e.g. AM < 5.0), the phenomenon 
that Jr decreases with K increasing does not occur. However, in case of(e.g. AM > 5.0), that 
phenomenon emerges and the maximum value of the average thermal current is gained at 
K = 1.0 for different mass gradients. It might be regarded as the result of the competition of 
the increment of thermal conductance originating from increasing [23] and the decrement 
of thermal conductance resulting from the existence of localized gradons modes|29] caused 
by mass gradient. In case of AM < 5, the rising thermal conductance induced by increment 
of K dominates the system. However, for increasing AM, localized gradons modes|29] is 


such non-ignorable that leads to the decrement of the thermal currents. 


Figure 7(a) depicts the average thermal currents (Jp, Jz) as a function of the system size 
N for K = 1.0, A= 1.0, Mmas = 10.0, Minin = 1.0, Ty = 0.5, Ty, = 0.1. It is found that 
the heat fluxes Jp/Jz, all decrease if the system size increases. Obviously, the difference of 
amplitude of J, and Jr gradually becomes smaller and smaller, which indicates thermal 
rectification effect gradually vanishes at a larger system size. 


According to above investigations on how the shape parameters influence the thermal 


FIG. 8: The heat fluxes (Jr, Jz) as a function of AT with Tr being fixed. The parameters are 
A =10.0, K = 1.0, Mmar = 20.0, Mmin = 1.0, Tr, = 0.1, N = 32. 


rectification efficiency, the optimum parameters are selected to construct the thermal recti- 


fier, whose rectification efficiency has reached R = 70 (see Fig. 8). 


IV. CONCLUDING REMARKS 


Heat conduction in the coupled rotator model with mass gradient has been investigation. 
Using nonequilibrium molecular dynamic simulation, the thermal currents through graded 
mass coupled rotator chain largely differ when the heat baths are reversed, which shows 
the thermal rectification effect. As for the coupled rotator model with graded mass, it 
is firstly found that the heat flux Jr in the direction of mass increasing is much larger 
than the heat current Jz in the direction of mass decreasing, which is opposite to that of 
the one-dimensional (1D) FPU lattices with graded mass[19]. It might be explained by 
match/mismatch of the frequency spectra. By further investigation, the larger is the mass 
gradient AM, the more evident is the thermal rectification effect, which is in accordance 
with the graded mass FPU-/ lattices[19, 20]. Secondly, the shape parameter A and K 
of nearest-neighbor interaction have influences on the thermal rectification. However, the 
phenomenon that the heat flux Jr decreases with K increasing is surprised. The reason is 
that the competition of the increasing thermal conductance arisen from the increment of K 
and localized gradons modes induced by uniform mass distribution. Thirdly, It is found that 


the thermal rectification efficiency becomes larger with increasing the average temperature. 
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Due to increment of the average temperature, the proportion of high-frequency phonons in 


entire phonons spectrum enhances. It leads to enhancement of the degree of mismatch of 


frequency spectrum at the two ends of atomic chain. Finally, based on above investigations, 


a relative optimum thermal rectifier has been constructed and its thermal rectification 


efficiency R =| Jr/ Jz |= 70. 
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